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Introductory Remarks 

The damped harmonic oscillator is arguably the simplest open quantum sys- 
tem worth studying. It is also of great practical importance because it is an 
essential ingredient in the theoretical description of many quantum-optical ex- 
periments. One can assume rather safely that the quantum master equation of 
the simple harmonic oscillator wouldn't be studied so extensively if it didn't 
play such a central role in the quantum theory of lasers and the masers. Not 
surprisingly, then, all major textbook accounts of theoretical quantum optics 
[§|,|,|j5}|.|,|.|jl^,|ll|,[|jl|Jl|,|l| contain a fair amount of detail about damped 
harmonic oscillators. Fock state representations or phase space functions of some 
sort are invariably employed in these treatments. 

The algebraic methods on which we'll focus here are quite different. They 
should be regarded as a supplement of, not as a replacement for, the traditional 
approaches. As always, every method has its advantages and its drawbacks: a 
particular problem can be technically demanding in one approach, but quite 
simple in another. This is, of course, also true for the algebraic method. We'll 
illustrate its technical power by a few typical examples for which the standard 
approaches would be quite awkward. 

1 First Lecture: Basics 

The evolution of a simple damped harmonic oscillator is governed by the master 
equation 

d 1 

— Qt = iLo[Qt,a)a] - ^'^{^ ^ {a^agt - 2aQta) + Qto) a) 

- -Au {aa^ Qt - 'ia) Qta + Qtaa)) , (1) 

where a\a are the ladder operators of the oscillator; uj is its natural (circular) 
frequency; A is the energy decay rate; and v is the number of thermal excitations 
in the steady state that the statistical operator Qt = Qt (a^ , a) approaches for very 
late times t. We'll have much to say about the properties of the solutions of (|l|), 
but first we'd like to give a physical derivation of this equation. 
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Fig. 1. A two-level atom traverses a high-quality cavity, coupling resonantly to a priv- 
ileged photon mode of the cavity. Prior to the interaction, there is some initial photon 
state in the resonator and the atom is either in the upper state of the pertinent transi- 
tion (on the left) or in the lower state (on the right). After the interaction, the transition 
degree of the atom and the photon degree of the cavity are entangled 



1.1 Physical Derivation of the Master Equation 

For this purpose we consider the following model. The oscillator is a mode of the 
quantized radiation field of an ideal resonator, so that excitations of this mode 
{vulgo photons) would stay in the resonator forever. In reality they have a finite 
lifetime, of course, and we describe this fact by letting the photons interact with 
atoms that pass through the resonator. As is depicted in Fig. |l|, these atoms are 
also of the simplest kind conceivable: they only have two levels which - another 
simplification - are separated in energy by fiuj^ the energy per photon in the 
privileged resonator mode. Incident atoms in the upper level (symbolically: :?:) 
will have a chance to deposit energy into the resonator, while those in the lower 
level (lii) will tend to remove energy. 

The evolution of the interacting atom-photon system is governed by the 
Hamilton operator 

H = hujo) a + huja'^ a — hg[aa) + cr^a) , (2) 

which goes with the name "resonant Jaynes-Cummings interaction in the rotat- 
ing-wave approximation" in the quantum-optical literature. It applies as long as 
the atom is inside the resonator and is replaced by 

-ffficc = fiLoa) a + hcja^a (3) 

before and after the period of interaction. Here and a are the atomic ladder 
operators, 

and g is the so-called Rabi frequency, the measure of the interaction strength. 
Note that a^a and ctct^ project to the upper and lower atomic states, 

respectively. 
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Fig. 2. An atom takes time r = L/v to traverse a cavity of lengtii L at speed ii 



The interaction term in 



7 = acr 

and iJfree is essentially the square of 7 since 



is a multiple of the coupling operator 

— n rr^ 



7^ = a' a + a'' a . 



So 



and 



can be rewritten as 

H = hcjj'^ — hgj 



free 



(6) 
(7) 
(8) 



That H commutes with -fffree and both are just simple functions of 7 will enable 
us to solve the equations of motion quite explicitly without much effort. 

We denote by Qt the statistical operator describing the combined atom-cavity 
system. It is a function of the dynamical variables a\ a, a\ a and has also a 
parametric dependence on t, indicated by the subscript, 



(9) 



Since a statistical operator has no total time dependence, Heisenberg's equation 
of motion, 

0=^,. = |p.4[p.i^], (10) 
becomes von Neumann's equation for the parametric time dependence, 

d 

Now suppose that t is the instant at which the atom enters the cavity; then it 
emerges at time t + r, and we have 



■[QuH] 



(11) 



--Ht 
= c ri gtcr^ 



e'^^gte-'^^ 



en 



(12) 



after we use (^ and introduce the abbreviation = gr. This phase is the 
accumulated Rabi angle and, for atoms moving classically through the cavity of 
length L with constant velocity v, we have t — L/v and — gL/v; see Fig. ^. 

Clearly, the [. . .] term in (|l^) accounts for the net effect of the interaction, 
that part of the evolution that happens in addition to the free evolution generated 
by -fffrcc- We have 



ji07 — cos (07) + i sin (^7) — cos ^^Vt^^ + 



sm 



(13) 
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which is just saying that the cosine is an even function and the sine is odd. 
Further we note that the identities 

F (72) = o-tcrF (aat) + (Tcr'fF (fll'a) , 

(72) = cTfltF (aat) + F (aa^) aa^ (14) 

hold for all functions ^"(7'^). They are immediate implications of familiar rela- 
tions such as af{a^a) = f{aa^)a, af{a''a) = cr/(l), and f{a''a) = 0-1^/(0). We 
use (^3|) and (Q) to arrive at 

gi^T = (tVcos (^i^VaaJ ^ + cra^ cos (^(j)VaJa^ 

sin ( (pV aat J sin ( 0-\/aaJ J 
+ icrat ^ ^ + i ^ ^ atjt . (15) 

In terms of the 2x2 matrix representation for the a's that is suggested in 
and (^), this has the compact form 



with the photon operators 

/ \ ~ / \ sin I (fiVaaJ J 

C = cos(<?!)\/aat) , C = cos ( 0\/at^ ) , S = ^ ^ , (17) 

^ ^ ^ ^ vaat 

and the adjoint of (|l^) reads 



C -iS-t 
-iS* C 



(18) 



We use these results for calculating the net effect of the interaction with one 
atom on the statistical operator gi(p'') of the photon state. Initially the total 
state Qt = Qt^^'^ 0t^^'^ is not entangled, it is a product of the statistical operators 

referring respectively to the photons by themselves and the atom by itself. At 

(ph 

t + T 



the final instant t + t, we get g[^^ by tracing over the two-level atom, 



gifj = Trat{ft+r} = e--™'«Trat {ei'^^giP'^^^,i^*)e~i'^^} e^™'" . (19) 

(at) 

To proceed further we need to specify the initial atomic state g). , and for the 
present purpose the two situations of Fig. |l| will do. 
On the left of Fig. |] we have atoms arriving, 

g[^'^ ^ |^)(^| = ata = ( J = ( J) (1,0) , (20) 
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and ( ]l9| ) tells us that 



t't+r 



Likewise, in the situation on the right-hand side of Fig. ^ we have 




1 



(0,1) 



(21) 



(22) 



and get 



(Ph) 
Qt+r = C 



—lujTa'a rp^ 



i^t 
c 



(23) 



We remember our goal of modeling the coupling of the photons to a reservoir, 
and therefore we want to identify the effect of very many atoms traversing the 
cavity (one by one) but with each atom coupled very weakly to the photons. 
Weak atom-photon interaction means a small value of so that only the terms 
of lowest order in (p will be relevant. Since the = version of both (|l|) and 
"231), that is: 



(ph) _ „—iuJTa^aApi^)A(JJTa''a 



Qi + T 



8t e 



(24) 



is just the free evolution of g'f^^\ the additional change in pj^'''' that results from 
a single atom is 



Qi 



(ph) 



, , , , sin ((f)V aflt ) sin ((j)y/aaJ ) 

UV^)gi''''^co8UV^) 

^ ^ ^ ^ vaa^ vaa^ 



9^' 



1 - 2 / t (ph) 



Jph) 



(25) 



for a atom arriving, and 

, , sin (0\/aat | , sinl^vaat ) 

ph) fi r~r\ , V / (ph t V / 

-^^ ' (0V aTfl 1 H — o,Qt " — ^ 



g{ cos I 



IS ^0\/ flta ^ 

(ataeiP'^) - 2ag^t^a^ + gi^'^^a^a) + 0(0*) 



(26) 
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for a zii atom. So, for weak atom-photon interaction the relevant terms in 
and ( p6| ) are the ones proportional to 0^. 

Atoms arriving at statistically independent times (Poissonian statistics for 
the arrival times) will thus induce a rate of change of g'f^^ that is given by 



d (ph) 

weak 

t (pli) r, t (ph) I (pli) t 

aa' Ql — 2a' Ql a + gl aa' 



(atae|P^) - 2ag'f''\^ + g^^'^^a^a) (27) 



where ri and r2 are the arrival rates for the atoms and the atoms, re- 
spectively. This is to say that during a period of duration T there will arrive on 
average riT atoms in state and atoms in state ij:. 

Since the weak interaction with many atoms is supposed to simulate the 
coupling to a thermal bath (temperature 6*) , these rates must be related to each 
other by a Maxwell-Boltzmann factor, 

where ^' > is a convenient parameterization of the temperature. Also for mat- 
ters of convenience, we introduce a rate parameter A by writing ri0^ — Av, 
r2(f>'^ = A{v -f 1) and arrive at 



d _ d 



d 



weak 



— iijj[gt, a] — -^^('^ + {'^^'^Qt — 2agta^ + gta^a) 

— 2^^^ {o-a^ Qt ~ 2a) gta + ptaa^) 

= Cgt , (29) 

where the replacement g^^^^ — > gt is made to simplify the notation from here 
on. It should be clear that the 0{(jy') of ( p5|) and (g^) terms are really negligible 
in the limiting situation of ri , r2 ^ A and cfP' <^ 1 with finite values for the 
products TicjP' and r20^. 

Equation (^) is, of course, the master equation (|l|) that we had wished to 
derive by some physical arguments or, at least, make plausible. From now on, 
we'll accept it as a candidate for describing a simple damped harmonic oscillator 
and study its implications. These implications as a whole are the ultimate justifi- 
cation for our conviction that very many crucial properties of damped oscillators 
are very well modeled by (p9|). 

Before turning to these implications, however, we should not fail to mention 
the obvious. The Liouville operator C of ( p9| ) is a linear operator: the identities 

C{\g) = \Cg , C{gi + g2) - Cgi + Cg2 (30) 

hold for all operators g^ gi, and g2 and all numbers A. 
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1.2 Some Simple Implications 

As a basic check of consistency let us first make sure that (^9|) is not in conflict 
with the normalization of Qt to unit total probability, that is: {gt} — \ for 
all t. Indeed, remembering the cyclic property of the trace, one easily verifies 
that 

^TT{Qt}^TT\^ot\^Tr{CQt} = Q, (31) 



At {dt 

as it should be. Much more difficult to answer is the question if (^9|) preserves 
the positivity of gt; we'll remark on that at the end of the third lecture (see 
Sect, m on p. H). 

Next, as a first application, we determine the time dependence of the expec- 
tation values of the ladder operator a\ a and the number operator a^a. Again, 
the cyclic property of the trace is the tool to be used, and we find 

- ^ Tr {atp,} ^ Tr {at|,,| . (i. _ lA){a^)^ , 

^(«M* = -^((«M*--)> (32) 
which are solved by 



(at), = (at)oe-^^/2ei-^ 
(a), = (a)oC-^^/2e-i-*, 
itfl)^ J, + ((ata)o - J^) e"^* , (33) 



respectively. Their long-time behavior, 

t^oo : (at)t ^ , {a)t ^ , (a^a)* ^ v , (34) 
seems to indicate that the evolution comes to a halt eventually. 



1.3 Steady State 

If this is indeed the case, then the master equation (^ ) must have a steady 
state g'^^'^\ As we see in (p2|), it is impossible for uj and A to compensate for 
each other and, therefore, must commute with the number operator a^a, 
and as a consequence is must be a function of a^a and cannot depend on and 
a individually. Upon writing /(a^a) for this function, we have 

= ^^g(^^^ = Cfia^'a) = -A{v + 1) [a^'afia'ia) - af{a'<a)a''] 

- Aiy[aa''f(a''a)-a''f{a''a)a] , (35) 
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and this implies the three-term recurrence relation 

{a^a + 1) [{i^ + l)f{a^a + 1) - t^/(ata)] = a^a [{i^ + l)/(a^a) - iyf{a^a - 1)] , 

(36) 

which, incidentally, is a statement about detailed balance (see for further 
details). In this equation, the left-hand side is obtained from the right-hand side 
by the replacement a^a ^ 0^0+ 1. Accordingly, the common value of both sides 
can be determined by evaluating the expression for any value that a^a may have, 
that is: for any of its eigenvalues (a'a)' = 0, 1, 2, . . . . We pick (a'a)' = and 
find that either side of (^6|) must vanish. The resulting two-term recursion, 

{,y + l)f{a^a)^i.f{a^a-l) (37) 

is immediately solved by f{a^a) — /(0)[z^/(z^-f 1)]°^^" and, after determining the 
value of /(O) by normalization, we arrive at 



This steady state is in fact a thermal state, as we see after re-introducing the 
temperature O of (p8|). 



1 — cxp — ■ 



e>^p(-#^«^«) • (39) 



Indeed, together with (jsj) this tells us that, as stated at (|l|), 'V is the number 
of thermal excitations in the steady state" . And the physical significance of A - 
it "is the energy decay rate" - is evident in (^). We might add that is the 
decay rate of the oscillator's amplitude (a), which is proportional to the strength 
of the electromagnetic field in the optical model. 

As the derivation shows, the steady state of ( [29| ) is unique, unless A — 0. 
Indeed, if A = but oj ^ 0, Cg ~ is solved by all g — f{a^a) irrespective of 
the actual form of the function /. We'll take A > for granted from here on. 

1.4 Action to the Left 



In (|32D we obtained differential equations for expectation values from the equa- 
tion of motion obeyed by the statistical operator, the master equation of (p9|). 
This can be done systematically. We begin with the expectation value of some 
observable X and its time derivative. 



{X),= Tr{Xg,} , ^(x)^ = Tr 



(40) 



and then use (p9|) to establish 



^^{X)^=Tr{XCgt} = {XC)^, (41) 
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where the last equation defines the meaning of XC, that is: the action of C to 
the left. The cyclic property of the trace is crucial once more in establishing the 
explicit expression 

XC = icj[ata, X] - ^A{i^ + 1) {Xa'fa - 2a'<Xa + a^aX) 

- ^Ai^ {Xaa'< - 2aXat + aa^ X) . (42) 



When applied to a\ a, and a^a, this reproduces (|3^), of course. 
How about (^l])? It is also contained, namely as the statement 

Tr {l£gt} = for aU gt, or 1£ = , (43) 

which is a statement about the identity operator 1. 

Homework Assignments 

1 Take the explicit forms of e^'t'^^ and e~^'^'>' in and verify that 

gi07 e-i-^T ^ 1 ^ e-'^^J e^'^T = 1 . (44) 



2 According to (Bq), the steady state (38) is a right eigenvector of C with 
eigenvalue 0, Cg^^^'' — 0. What is the corresponding left eigenvector g^^^'^ 
such that ^("")£ = 0? 

3 Reconsider (|3|) and (^). Show that these equations identify some other 
eigenvalues of £ and their left eigenvectors. 

4 Use the ansatz 

gt = m[l--\it)f'' (45) 

in the master equation (|29|), where X{t) X = 1/(1 + z/) for t ^ oo. 
Derive a differential equation for the numerical function X{t), and solve it 
for arbitrary A(0). [If necessary, impose restrictions on A(0).] Then recognize 
that the solution reveals to you some eigenvalues of C. Optional: Identify the 
corresponding right eigenvectors of C. 

2 Second Lecture: Eigenvalues and Eigenvectors of C 
2.1 A Simple Case First 

When taking care of homework assignment ^, the reader used (^) to establish 
Cg,^-A^^-^^^[Xa^a-{l-X)]gt (46) 
for gt of ( p5| ) and found 
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by differentiation. Accordingly, (Esl) solves (p^) if X{t) obeys 



1 dA _ d 1 
A2 di ^ di A 



which is solved by 



m = 



(j. + 1) - [(^. + 1) - 1/A(0)] c 



(48) 



(49) 



where the restriction A(0) > is sufficient to avoid ill-defined values of A(t) at 
later times, and A(0) < 1 ensures a positive gt throughout. 
With (E|) in (||) we have 



9t 



E 

n=0 



p-nAt JO) 



(50) 



where the gl^'^^s are some functions of a^a. In particular, ~ g'^^^'' is the 
steady state ( p8| ) that is reached ior t ^ oo when X{t) — > + 1). Since ( pO| ) 
is a solution of the master equation (p9|) by construction, it follows that 



(51) 



n=0 



n=0 



holds for alH > and, therefore, pi"'' is a right eigenvector of C with eigenvalue 
-nA, 

= . (52) 

As defined in ( |5C| ) with of ( ^ ) and X{t) of ( |49| ) on the left-hand side, the gn ^'s 
depend on the particular value for A(0), a dependence of no relevance. We get 
rid of it by introducing a more appropriate expansion parameter x in accordance 
with 



1 



m 



{v + 1)(1 + x) or X 



1 



(j. + l)A(0) 



1 



-At 



(53) 



so that counting powers of e is done by counting powers of x. Then 
1 



{v + l){l + x) 



1 - 



H \ a a oo 

1){1 + X) 2^^ 



(54) 



is a generating function for the g\i^^s with the spurious A(0) dependence removed. 

The left-hand side of ( 54 ) can be expanded in powers of x for any value of 
V > but we'll be content with a look at the v = Q case and use a different 
method in Sect. 2.2 to handle the general situation. For v = the power series 
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(p4[) is simplicity itself, 



\ X 



0: 5:x>(") = 



1 + a; V 1 + a; 

n=0 



m=0 ^ ' 

f;(-ir("'"+™)X«+-, (55) 



so that 



. = 0: = M, (56) 

It is a matter of inspection to verify that 



: (57) 

are the v — Q right eigenvectors of C to eigenvalues 0, —A, — 2A, .... 
We obtain the corresponding left eigenvectors from 

oo 

^. = 0: (l + y)"'"-^^"^™ (58) 
after verifying that this is a generating function indeed. For v = Q, (|4|) says 

(1 + yf^^C ^ -Aa^a{l + y)"^" + Aa\l + y)"^"a 
= -Ayata(l+yr'"-l 

= -Ay^(l + y)"'« (59) 
and the eigenvector equation 

e'^^^C = -mApt^) (60) 

gives 

oo oo oo 

Y: v'-st^C = V-i-mAYgt^ = ^Ay^ ^ y^g^^ , (61) 

m=0 m=0 ^ m=0 



and now (^ and (|6l| ) establish (|5q). So we find 

. = 0: (62) 
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of which the first few are 

-(0) t 

g\ ~ a' a , 

4°) = lata(ata- 1) . (63) 

For n = and n = 1 this just repeats what was learned in ( ^3| ) and (^) (recall 
homework assignments 2 and 3). 

As dual eigenvector sets, the Qn'^'s and ^m^'s must be orthogonal if n 7^ m, 
which is here a statement about the trace of their product. It is simplest to deal 
with them as sets, and we use the generating functions to establish 



Y: 2/'"Tr{^(°),(°)}x"= Tr (l + 2;r'«^(- 

■n,n=0 [ ^ 

l + x ^ V ^'l + x 



1 ~ xy ^ — ' 



= y^-^m.na;", (64) 



from which 



Tr 



{q^} Q^^^] = 5ra,n (65) 

follows immediately. This states the orthogonality of the v = Q eigenvectors and 
also reveals the sense in which we have normalized them. 

In the third lecture we'll convince ourselves of the completeness of the eigen- 
vector sets. Let us take this later insight for granted. Then we can write any 
given initial statistical operator pt=o = /(a^i) as a sum of the Qn\ 

oo 

e*=o = E"»^^'n^ (66) 



n=0 



and solve the master equation (G3) by 



^are-^^p(«). (67) 



Qt = 

11=0 

As a consequence of the orthogonality relation (^5|), we get the coefficients ai"'' 
as 

Tr{ft(°)g,=o} . (68) 
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For V = 0, in particular, they are 



which tells us that (|66| ) and (|6^) are expansions in moments of the number 
operator. Put differently, the identity 

00 

/(ata) = ^gWTr{ft(°)/(ata)} (70) 

n=0 

holds for any function f{a^a) that has finite moments. For most of the others, 
one can exchange the roles of Qn^ and gl^"^ , 

00 

/(ata) = 5]Tr{/(ata)^?(")}^i°). (71) 

71=0 

A useful rule of thumb is to employ expansion (^0|) for functions that have the 
basic characteristics of statistical operators [the traces in ( |70| ) are finite], and 
use ( |7l| ) if /(a^a) is of the kind that is typical for observables [such as a^a for 
which the n = trace in ^u) is infinite, for example]. 



2.2 The General Case 



Let us observe that some of the expressions in Sect. 2.1 are of a somewhat simpler 
structure when written in normally ordered form - all operators to the left of 
all a operators - as exemplified by 

A(l- A)"^" = A : e-^«^" :, 
(l + y)»^« = : e?/"^" : , 



p = : - ] = —: ■■ (a^a) : = ^ • (72) 

m J m! m! 



We regard this as an invitation to generalize the ansatz ( [45| ) to 

g^ = :—e~ l-"^ ~ - " W] : (73) 

K{t) 



where the switching from A to k = 1/A is strongly suggested by (48). Note 
that for the following it is not required that a*{t) is the complex conjugate of 
a{t), it is more systematic to regard them as independent variables. We pay due 
attention to the ordering and obtain 
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for the parametric time derivative of Qt- The evaluation of Cgt is equahy straight- 
forward once we note that a^'s on the right and a's on the left of Qt are moved 
to their natural side with the aid of these rules: 

d 1 

Qto^ = Qt + [Qt,a)] = a) Qt + TrBt = a^Pt (a^ - ce*)Qt , 

oa K 

d 1 

aQt = Qta + [a, Qt] = fiija + = fta - ^ Q?) C''^) 

of which ^ ^ 

[gt, a^a] = [a) ~ a*)Qta H — a^f?t(a — a) (76) 

K K 

is an immediate application. Upon equating the numerical coefficients of both 
Qt and (a^ ^ a*)Qt{a — a) we then get a single equation for 

^ = + , (77) 

and the coefficients of (a^ — a*)£it and Qt{a — a) supply equations for a{t) and 
a*(i), 

We have, in fact, met these differential equations before, namely (^7|) in ( p8| ) and 
in (^). Their solutions 

K{t) = 1/ + 1 - (j/ + 1 - Ko) e""^* , 
a{t) = aoe-i^^e-^*/2 , 
a*(t) =aSei'^*e-^^/2 , (79) 

where kq, ao, are the arbitrary initial values at i = [not to be confused 
with the time-independent coefficients an^ of (|66|)-(|69|)], tell us that the time 
dependence of 

Qt = e^^Qt^o (80) 

contains powers of e~^* combined with powers of e^^'^* ~ Therefore, the 

values 

=-ikuj- {n+^\k\)A with n = 0,1,2,.. . and fc = 0,±1, ±2,.. . (81) 

must be among the eigenvalues of C. In fact, these are all eigenvalues, and the 
expansion of the generating function ( |73| ) 



E E 6We"^^^*-(^+^1^0^*^,W (82) 



OO OO 

Qt = 

n—0 k— — oo 



yields all right eigenvectors of £, 

jOgi'^ = [-ikco-{n+^\k\)A]Qi^K (83) 
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We'll justify the assertion that these are all in the third lecture. Right now we 

(k) 

just report the explicit expressions that one obtains for the eigenvectors g„ and 
the coefficients bi^^ . They are 0@ 

(4_„ti(l^l+fc) :L(l^l)f^)e"^: ~ k) (34) 

and 

{n + \k\)l \iy + 1 J "° " \j, + l~KoJ 

(85) 

where the L„ s are Laguerre polynomials. Note that all memory about the 
initial values of k, a, and a* is stored in bl^\ as it should be. The right eigenvec- 
tors of (^) and the left eigenvectors of (^) constitute the two so-called 
damping bases associated with the Liouville operator C of (^9|). 



Homework Assignments 

5 Verify that (pih reproduces (poj for fc = and ly ^ 0. 



6 Show that (|82|) with (|8^) and (|5|) is correct. What is b^n^ for = + 1, 
ao = {v + I)f3, and ajj = (i^ + l)/3*? [You need some familiarity with Bessel 
functions, Laguerre polynomials, and relations between them; consult the 
Appendix if necessary.] 

7 Consider 

U = e""^ e"*" e^ = e"*'^ e""^ e^ " = e""^ + e^ " (86) 

with a*{t), a{t), X{t) such that dU/dt = UC. Find the differential equations 
obeyed by a*, a, A and solve them. 

8 Show that this also establishes the eigenvalues ( |8l| ) of C. 

(k) 

9 Extract the left eigenvectors p„ of C. Normalize them such that 

^" [l + iyj (n+|fc|)! ■''^ { ly J ■ ■ 

What do you get in the limit ly ^ 07 

10 State g'i!'^ explicitly for fc = 0, 71 = 0, 1,2 and for fc = ±1,±2, n = 0. 
Compare with (|63|). 

11 Use the two generating functions (^) and (p6h to demonstrate 
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3 Third Lecture: Completeness of the Damping Bases 

3.1 Phase Space Functions 

As a preparation we first consider the standard phase space functions f{Q', P') 
associated with an operator F{Q,P) where Q,P are a Heisenberg pair (col- 
loquially: position Q and momentum P). This is to say that they obey the 
commutation relation 

[Q,P]=i (89) 
and have complete, orthonormal sets of eigenkets and eigenbras, 

Q\Q') = \Q')Q' , P\P') = \P')P' , 

{Q'\Q = Q'iQ'l , {P'\P = P'{P'\ , 

{Q'\Q") = S{Q' - Q") , {P'\P") = d{P' - P") , 

\Q'){Q'\ = S{Q - Q') , \P'){P'\ = S{P - P') , 

J dQ' \Q'){Q'\ = 1 , JdP' \P'){P'\ = 1 , (90) 

where Q',Q" and P',P" denote eigenvalues. The familiar plane waves 

iQ'in = ^ (91) 



relate the eigenvector sets to each other. 

By using both completeness relations, we can put any F = F{P, Q) into its 
Q, P-ordered form. 



P(P,Q) = JdQ'dP'\Q'){Q'\F\P'){P'\ 

■{Q'\F\P') 



dQ'dP'\Q'){Q'\ 



\P'){P'\ 



= jdQ'dP'5{Q-Q')f{Q',P')5{P-P') 



= f^Q'PK. . ,^m-^p), (92) 

<y,F— ordered 

where the last step defines the meaning of the semicolon in f{Q; P). Thus, the 
procedure is this: evaluate the normalized matrix element 

;W^n=<«Hn=TJf^l, (93) 



(Q'|P') I {Q'\P') 

then replace Q' ^ Q, P' ^ P with due attention to their order in products 
- all Q's must stand to the left of all P's - and so obtain F = f{Q;P), the 
Q, P-ordered form of F. The P, Q-ordered version of F is found by an analogous 
procedure with the roles of Q and P interchanged. 
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The fraction that appears in the trace formula of (93) is equal to its square, 



\P'){Q'\ 



\P'){Q'\ 

mp') 



(94) 



but, not being Hermitian, it is not a projector. It has, however, much in common 
with projectors, and this is emphasized by using 



to turn it into 



= 27:\P'){P'\Q'){Q'\ ^ 2^5{P - P')5{Q - Q') 

which is essentially the product of two projectors. Then, 

/(Q', P') = Tr {F 27t6{P ~ P')5{Q ~ Q')} 

is yet another way of presenting /(Q', P'). 

In ( p^ ) and (|9^) we recognize two basis sets of operators, 

B{Q', P') = 2^6{Q - Q')6{P - P') , 
B{Q', P') = 2^5{P - P')S{Q - Q') , 



(95) 



(96) 



(97) 



(98) 



labeled by the phase space variables Q' and P' . Their dual roles are exhibited 
by the compositions of (^2|) and (p7|). 



F{Q,P) = J^^^B{Q',P')Tt{fB{Q',P')} 



and 



f{Q',P')= Tr\B{Q',P' 



dQ" dP" 
2^ 



f{Q",P")BiQ",P") 



and therefore 



Tr {b{Q', P')B{Q'\ P")} = 2TTSiQ' - Q")S{P' - P") 

states both their orthogonality and their completeness. 
We note that the displacements 



Q^Q-Q' 



P ^ P- P' 



(99) 



(100) 



(101) 



(102) 



that map out all of phase space (see Fig. ||) are unitary operations, and so is the 
interchange 

Q^P, P^-Q. (103) 



18 



Berthold-Georg Englert and Giovanna Morigi 



P' 



•{Q',P') 

>g' 



\seed S(0, 0) 

Fig. 3. The unitary displacements Q-^Q — Q',P^P~P' map out all ol phase 
space, turning the basis seed B{0,0) into all other operators B{Q',P') of the basis 



As a consequence, B{Q\P') and B{Q',P') are unitarily equivalent to their re- 
spective basis seeds 

B(0,0) = 2T:d{Q)d{P) and 5(0,0) = 27r6{P)d{Q) , (104) 

and these seeds themselves are unitarily equivalent and are also adjoints of each 
other, 

B(0,0)t = 5(0,0) , 5(0,0)^=5(0,0), (105) 

Here, 5(0, 0) is stated as a Q, P-ordered operator and 5(0, 0) is P, Q-ordered. 
The reverse orderings are also available, as we illustrate for 5(0,0): 



(P'|5(0,0)|Q') ^ ^JP'\Q"){Q"\P"){P"\Q') 



{P'\Q') {P'\Q' 

giving 

OO 

kl 



= e'PQ , (106) 

'=0, P"=0 



OO . ^ 

5(0, 0) = c^^' Q = ^ iyP'^Q^ (107) 

fc=0 



where we meet a typical ordered exponential operator function. We get the Q, P- 
ordered version of 5(0, 0), 

5(0,0) = e-i<5;-P ^ (108) 

by taking the adjoint. 

These observations can be generalized in a simple and straightforward man- 
ner [^. The operators 

A ei^^' Q = \ e-i^Q; P , A e-^^Q; = A e^^^^ Q 

with A, A > 1 and AA = A + A (109) 

are seeds of a dual pair of bases for each choice of A, A because the orthogonality- 
completeness relation 

Tr { A ei^(^ - ' - A c-i^(« - ' - } = 2^ W - g")<5(P' - P") 

(110) 
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holds generally, not just for A = 1, A = oo as we've seen in (101 ). Wc demonstrate 
the case by using the two completeness relations of ( pO[ ) to turn the operators into 
numbers, and then recognize two Fourier representation of Dirac's 6 function: 



Tr 



{■ I 



2tt 



Xe 



-iA(Q 



'){P-P") 



J 2Tr J 2tt 

= 2TTe'^^Q'P' - Q"P"h{Q' - Q")6{P' - P") 
= 2it5{Q' - Q")5{P' - P") , 



(111) 



indeed. The equivalence of the A and A versions of (109) is the subject matter 
of homework assignment |l^. 

The basis seeds (109) can be characterized by the similarity transformations 
they generate, 

Q X(^^P;Q ^\(^^P\Q {\-X)Q or Q^(1-A)Q, 

P Xe^^P'^Q = Xe^^^'^Q (l-X)P or P ^ (1 - A)P , (112) 

which are scaling transformations essentially (but not quite because 1 — A = 



1/(1 — A) < and the cases A = 1 or A 
with the aid of identities such as 

1 d 



1 are particular). One verifies (112) 



1 oQ 



P e 



i\P-Q 



(113) 



Equations (112) by themselves determine the seed only up to an over-all factor, 
and this ambiguity is removed by imposing the normalization to unit trace. 



Tr 



'^Q_^y^^iXPQ 



27r 



1 



(114) 



The cases A = 1, A — > oo and A — > oo, A = 1 are just the bases associated 
with the Q, P-ordered and P, Q-ordered phase space functions that we discussed 
above. Of particular interest is also the symmetric case of A = A = 2 which has 
the unique property that the two bases are really just one: The operator basis 
underlying Wigner's phase space function. Its seed (not seeds!) is Hermitian, 
since A = A — 2 in ( |l09| ) implies 



2gi2P;Q ^ 2e-^2Q;P ^ Uq^'^P'^Q 



t 



(115) 



It then follows, for example, that F — F'^ has a real Wigner function; in fact, 
all of the well known properties of the much studied Wigner functions can be 
derived rather directly from the properties of this seed (see for details). 

For our immediate purpose we just need to know the following. When A = 
A = 2, the transformation ( 112 ) is the inversion 



P-> -P 



(116) 
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For at = 2 ^/^(Q — iP), a = 2 ^/^{Q + iP) this means that the number oper- 



ator a^a — ■^{Q^ + — 1) is invariant. Put differently, the Wigner seed (115) 



commutes with a^a, it is a function of a^a: 2e^^^''*5 — f{o)a). For a, the in- 
version ( |116| ) requires af{a^a) — —f(a^a)a, and this combines with af{a^a) = 
f{a^a + l)a to teU us that f{a^a + 1) = — /(a^a). We conchide that 



2gi2P;Q ^ 2(-l)"^" = : 2e~2ata . 



after using the normahzation (114) to determine the prefactor of 2, 

2 

n=0 n=0 l>x~*l 



OO C30 

Tr {2(-l)"^'^| = 2 ^(-1)" = 2 ^(-x)" 



1 + a; 



(117) 



(118) 



The last, normally ordered, version in (117) of the Wigner seed shows the A = 2 
case of (^). 

Perhaps the first to note the intimate connection between the Wigner func- 
tion and the inversion ( |116| ) and thus to recognize that the Wigner seed is (twice) 

the parity operator (—1)°" " was Royer In the equivalent language of the 
Weyl quantization scheme the analogous observation was made a bit earlier by 
Grossmann A systematic study from the viewpoint of operator bases is 
given in Q . 



We use the latter form (117) to write an operator F{a\a) in terms of its 
Wigner function f{z*, z), 



F{a\a 



t n\ - 



dQ' dP' 
2^ 



f{z*,z) :2e 



-2 flf 



(119) 



where z* = 2-'^^^(Q' - iP') and z = 2-'^^'^{Q' + iP') are understood, and 

f{z*,z) ^ Ti I F{a\ a) : 2 e'^K " (° " ^) A (120) 



reminds us of how we get the phase space function by tracing the products with 
the operators of the dual basis (which, we repeat, is identical to the expansion 
basis in the A = A = 2 case of the Wigner basis). 



3.2 Completeness of the Eigenvectors of C 

The stage is now set for a demonstration of the completeness of the damping 



bases of Sect. |2^. We'll deal explicitly with the right eigenvectors gn of the 
Liouville operator £ of (|29| ) as obtained in ( |8^ ) by expanding the generating 
function ( |7^ ) . Consider some arbitrary initial state gt=o and its Wigner function 
representation 

f dQ'dP' ^ , ^ -2(a^ -z*)(a~z) x 

Qt=o= — ^ — p(z,z):2e V A J-., (121) 
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where p{z*,z) is the Wigner phase space function of Qt=o- Then, according to 
Sect. |2.2|, we have at any later time 



J 27r K.{t) 

with 

Ko = I , ao = 2: , tto = z* (123) 
in (|7|). In conjunction with (||), (||), and (||) this gives 

^* = E E /3We-^^^*-("+2lfc|)^*^,W (124) 

n— fc— — 00 



with 



-\k\)\ 

dQ' dP' 
"hi 



(125) 



But this is just to say that any given Qt has an expansion in terms of the gii 
for all t - any Qt can be expanded in the right damping basis. Equation ( |124[ ) 
also confirms that the exponentials exp(— ifccjt — (n + i|fc|)j4t) = exp(A^'^ 
constitute all possible time dependences that Qt might have. Clearly, then, the 
A^'^^'s of ( ^ ) are all eigenvalues of £, indeed, and the right eigenvectors Qn'' of 
(HI) are complete. The completeness of the left eigenvectors Q'n'^ of ( ^7| ) can be 
shown similarly, or can be inferred from (p8|). 

Note that this argument does not use any of the particular properties that 
Qt=Q might have as a statistical operator. All that is really required is that 
its Wigner function p[z*,z) exists, and this is almost no requirement at all, 
because only operators that are singularly pathological may not possess a Wigner 
function. Such exceptions are of no interest to the physicist. 

More critical are those operators Qt=o for which the expansion is of a more 

(k) 

formal character because the resulting coefficients fin are distributions, rather 
than numerical functions, of the parameters that are implicit in p{z*^z). In 
this situation the recommended procedure is to expand in terms of the left 
eigenvectors Qii instead of the right eigenvectors Qn ■ 

The demonstration of completeness given here relies on machinery developed 



in 17, 18| , |l9| . An alternative approach can be found in [g^, where the case w = 0, 
1/ = is treated, but it should be possible to use the method for the general case 
as well. 



Note also that the time dependence in (122) is solely carried by the operator 
basis, not by the phase space function. This is reminiscent of - and in fact closely 
related to - the interaction-picture formalism of unitary quantum evolutions. 
Owing to the non- unitary terms in £, those proportional to the decay rate A, 
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the evolution of gt is not unitary in (122). Of course, one could also have a 



description in which the operator basis does not change in time, but the phase 
space function does. It then obeys a partial differential equation of the Fokker- 
Planck type. Concerning these matters, the reader should consult the standard 
quantum optics literature as well as special focus books such as p^ . 

3.3 Positivity Conservation 

Let us now return to the question that we left in limbo in Sect. |l.2| : Does the 
master equation (p9| ) preserve the positivity of gt? 

Suppose that gt=o is not a general operator but really a statistical operator. 



Then the coefficients Pn'^ in (124) are such that the right-hand side is non- 
negative for i = and has unit trace. Since Pq"' is the steady state g^-^^^ of ( p8| ) 
and ^q"' = 1 is the identity, we have 

Tr { gl'^ } = Tr { g^^^ ) } - <5„,o4,o (126) 

as a consequence of (|8|), and 

/3« = Tr{ftWft.o} (127) 



implies /Sg"'' = Tt {gt^o} — 1- So, all terms in (124) are traceless with the sole 



exception of the time independent n = 0, fc = term, which has unit trace. This 



demonstrates once more that the trace of gt is conserved, as we noted in Sect. |l.2 
already. 



The time independent n = 0, fc = term is also the only one in (124) that 



is non-negative by itself. By contrast, the expectation values of git'' are both 
positive and negative if n > 0, fc = and even complex if fc 7^ 0. Now, since 
fft=o > 0, the n = 0,k = term clearly dominates all others for t = 0, and then it 
dominates them even more for i > because the weight of pg°'' = g^^^'^ is constant 

(k) 

in time while the other gn 's have weight factors that decrease exponentially 



with t. We can therefore safely infer that the master equation (^9|) conserves the 
positivity of gt ■ 

3.4 Lindblad Form of Liouville Operators 

The Liouville operator of (|2^) can be rewritten as 

Cg — iLLi[g, fltfl] -I- + 1) ([fl' 0^^] + ['^Qi "■^]) 

+ ^Aiy{[a^,ga] + [a^g,a]) (128) 

and then it is an example of the so-called Lindblad form of Liouville operators. 
The general Lindblad form is 
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where H = H'^ is the Hermitian Hamihon operator that generates the unitary 
part of the evolution. Clearly, any C of this form will conserve the trace, but the 
requirement of trace conservation would also be met if the sum were subtracted 



in (129) rather than added. As Lindblad showed however, this option is only 
apparent: all terms of the form , gVj] + [Vj g, Vj] in a Liouville operator must 
come with a positive weight. He further demonstrated that all £'s of the form 
(129) surely conserve the positivity of g provided that all the 1^ 's are bounded. 



This fact has become known as the Lindblad theorem. In the case that some Vj 
is not bounded, positivity may be conserved or not. A proof of the Lindblad 
theorem is far beyond the scope of these lectures. Reference is perhaps a 
good starting point for the reader who wishes to learn more about these matters. 
We must in fact recognize that (|128D obtains for 



Vi ^ ^ lA{iy + 1) a'' , V2^^J\Ava (130) 

in ( |129| ) and these are actually not bounded so that the Lindblad theorem does 
not apply. Fortunately, we have other arguments at hand, namely the ones of 
Sect. |3.3| , to show that the master equation ( p9| ) conserves positivity. 

But while we are at it, let's just give a little demonstration of what goes 
wrong when a master equation is not of the Lindblad form. Take 

^Ot = VV^gt~2V^gtV + gtVV^ , (131) 
at 

for example, and assume that the initial state is pure, gt=o — |V'o)(V'o|- Then we 
have at t = di > 



dgt 

Qt=dt = Qt=o+Cit— 
ot 



|Vo>(^o|+dt(|V2)(Vo|-2|Vi)(V'i| + IV'o)(^2|) , (132) 



t=0 



where = l^^lV'o) and 1-02) = T^F'''|?/'o). Choose IV'o) such that (V'i|^/'o) = 
and (V'ilV'2) = 0, and calculate the probability for at time t = dt, 



(V'l|£'t=dt|'0 



(V'llV'i) 



— = -2dt < ! (133) 



Positivity is violated already after an infinitesimal time step and, therefore, (131) 
is not really a master equation. 



Homework Assignments 

12 Show the equivalence of the A and A versions of ( |109| ). 

13 If F = fi{Q;P) = f2{P;Q), what is the relation between /i(Q',P') and 
/2(P',Q')? 

14 Consider F — ^{QP + PQ), the generator of scaling transformations. For 

real numbers /z, write e^^^ in Q, P-ordered and P, Q-ordered form. 

15 Find the Wigner function of the number operator a^a. 

16 Construct an exphcit (and simple!) example for (131)-(133), that is: specify 
V and \ipo)- 
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4 Fourth Lecture: Quantum- Optical Applications 

Single atoms are routinely passed through high-quality resonators in experiments 
performed in Garching and Paris, much like it is depicted in Fig. |l|; see the review 
articles [^,^,^ 29 3^ and the references cited in them. In a set-up typical of 
the Garching experiments, the atoms deposit energy into the resonator and so 
compensate for the losses that result from dissipation. The intervals between 
the atoms are usually so large that an atom is long gone before the next one 
comes, so that at any time at most one atom is inside the cavity, all other 
cases being extremely rare. And, therefore, the fitting name "one-atom maser" 
or "micromaser" has been coined for this system. 

The properties of the steady state of the radiation field that is established in 
the one-atom maser are determined by the values of several parameters of which 
the photon decay rate A, the thermal photon number v, and the atomic arrival 
rate r are the most important ones. Rare exceptions aside, an atom is entangled 
with the photon mode after emerging from the cavity and it becomes entangled 
with the next atom after that has traversed the resonator. As a consequence, 
measurements on the exiting atoms reveal intriguing correlations which are the 
primary source of information of the photon state inside the cavity. A wealth of 
phenomena has been studied in these experiments over the last 10-15 years, and 
wc look forward to seeing many more exciting results in the future. 



4.1 Periodically Driven Damped Oscillator 

To get an idea of the theoretical description of such experiments and to show 
a simple, yet typical, application of the damping bases of Sect. || without, how- 
ever, getting into too much realistic detail, we'll consider the following some- 
what idealized scenario. The atoms come at the regularly spaced instants t = 
. . . , — 2r, — T, 0, T, 2T, . . . and each atom effects a quasi-instantaneous change of 
the field of the cavity mode (= the oscillator) that is specified by a "kick opera- 
tor" /C which is such that 



QjT^o QjT+0 = (1 + K,)QjT-a (134) 

states how Qt changes abruptly a,t t = jT [j — 0,±1,±2, . . .). Physically speak- 
ing, this abruptness just means that the time spent by an atom inside the res- 
onator is very short on the scale set by the decay constant A. 

The evolution between the kicks follows the master equation (p9|). We have 

d 

— Qt^Cgt+IC^S{t-jT)gt-o (135) 



as a formal master equation that incorporates the kicks (134). The Cg^ term 
accounts for the decay of the photon field in the resonator. In Sect. 1.1 we 
derived the form of C by pretending that the dissipation is the result of a very 
weak interaction with very many atoms. But, of course, this model must not be 
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1 2 ^ At ^ ^ 

Fig. 4. Mean number of excitations of a periodically kicked oscillator; see text 



regarded as the true physical mechanism. Actually, the loss of electromagnetic 
energy is partly caused by leakage through the cavity openings (through which 
the atoms enter and leave) and partly by the ohmic resistance from which the 
currents suffer that are induced on the surface of the conducting walls. The 
ohmic losses are kept very small by fabricating the cavity from superconducting 
metal (niobium below its critical temperature). 

Thus, the term Cgt in ( |l35| ) has nothing to do with the atoms that the 
experimenter passes through the resonator in a micromaser experiment. These 
atoms interact strongly with the photons and give rise to the kick oper ator /C. 
Figure ^ shows what to expect under the circumstances to which ( |l35| ) refers. 

The mean number of excitations, (a^a)^, decays between the kicks (solid 

line) in accordance with (^) and changes abruptly when a kick happens (vertical 

dashed lines at At = 0.4, 0.8, 1.2, . . .). After an initial period, which lasts 

about a dozen kicks in Fig. ^, the cyclically steady sta te g ^^^^^^ is reached whose 



defining property is that it is the periodic solution of (135) 



(ess) 
St+T 



Bi 



(ess) 



Its value just before a kick is determined by 



Bt 



(ess) 



'o=e^^,tto=e^^(l+^)e„ 



and we have 



(ess) Ct (ess) 

Bt - ^ bU+o = c 



-C{T-t)(css) 

Bt=-o 



(136) 



(137) 



(138) 



for < t < r, that is: between two successive kicks. 



26 



Berthold-Georg Englert and Giovanna Morigi 



In passing, it is worth noting that a recent experiment [ pT| , in which photon 
states of a definite photon number (Fock states) were prepared in a micromaser, 
used a periodic scheme for pumping and probing. The theoretical analysis |Q 
benefitted from damping-bases techniques. 

The fine detail that we see in Fig. |^ is usually not of primary interest, partly 
because experiments tend to not resolve it. For example, if one asks how long 
it takes to reach the cyclically steady state, all one needs to know is the time- 
averaged behavior of the smooth dash-dotted — • — • — lines in Fig. ^ In the 
cyclically steady state, the meaning of "time-averaged" is hardly ambiguous, we 
simply have 

-gi-)^l£dtg[-\ (139) 

where it does not matter over which time interval we average as long as it covers 
one or more periods of g'f'^^^. The periodicity ( |l36D of g'f^^^ imphes that dgl^^^ /dt 
is zero on average and, therefore, we obtain 

+ \^JCQt% = (140) 



when time averaging ( |l35| ). When combined with what we get upon using (138) 
in B), 



1 - e-^r 



it yields the equation that determines p'^^*'-' , 



+ IC- ^TcT ^^"^ = . (142) 



1 - e- 



For /C = 0, it is of course solved by p^'^^^ = p^Q^ of (|3^). 

A master equation for the time-averaged evolution, that is: an equation 
obeyed by the time-averaged statistical operator gt, cannot be derived from 



(135) for the same reasons for which one cannot derive the macroscopic Maxwell 
equations from the microscopic ones. But they can be inferred with physical 
arguments that are more than just reasonably convincing. The task is actually 
easier here because we have to deal with temporal averages only whereas one 
also needs spatial averages in the case of electromagnetism. 
Imagine, then, that a linear time average is taken of (|135|), 



^^gt=Cgt+ICi?)gt , (143) 



where (?) gt is the ill-d eterm ined average of the summation in ( 135 ) that accounts 
for the periodic kicksJ134). We require, of course, that g'-^^'' is the steady state 



of (|143| ). In view of ( |142|) , this requirement settles the issue [ p3| . 
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which we now accept as the master equation that describes the time-averaged 
evolution. This is another case where an equation is ultimately justified by its 
consequences. 

We run a simple, but important consistency check on (144). If the spacing 
T between the atoms decreases, T — 0, and also the effect of a single atom, 
/C = pA4 with p —> 0^ such that their ratio r = p/T is constant, then the 
situation should be equivalent to that of Poissonian arrival statistics with rate 
r and each atom effecting a kick A4. Indeed, ( |l44|) turns into 

d _ 

— Qt^Cgt+rMgt, (145) 

as it should, because this is the familiar Scully-Lamb equation that is kn own 
to apply in the case of Poissonian statistics. Thus, with /C = rTM in (144) we 
obtain a master equation, 

d _ _ CT 

— gt = Cgt+rM ^t^P* , (146) 

at 1 — c '-^ 

that interpolates between the Poissonian Scully-Lamb limit of T = and that 

of highly regular arrival times, T = 1/r. 

The damping bases associated with C are the crucial tool for handling (144). 
We write 

oo oo 

B^ = Y.Y. (147) 

n— k— — oo 

— (k) 

and obtain differential equations for the numerical coefficients aii , 

aW(i)= Tr{^('='ft} , (148) 

by exploiting 

f{C)gi!:^ = /(aW)^,W , = -iku; - (n + ^\k\)A , (149) 

which holds for any function /(£) simply because gi^"^ is the right eigenvector 
of C to eigenvalue aI'^'' . In this way, (|144|) implies 



n',k' 1 - e~ ^ 



where 



{e^icgi':'} (151) 

is the matrix representation of the kick operator JC in the damping bases. The 
seemingly troublesome ratio in (144) is not a big deal anymore in (150), where 
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of course. 

Let us illustrate this for the particularly simple kick operator specified by 
JCg—p(a^ :g g — g ) with 0<p<l, (153) 



which describes the over-idealized situation in which an atom adds one photon 
with probability p and does nothing with probability 1 — p. Here, 

^ix'^ = Sk.k'ici^'} (154) 

so that the evolution does not mix s of different k values. As a further 
simplification it is therefore permissible to just consider the fc = terms. For 
1/ = (homework assignment |l^ deals with > 0) , the generating functions ( p5|) 
and ( pSi ) give 



m,n— 

1 — xv ^ — ' 

^ n=0 

where we let JC act to the left, 

{l + yf-IC^py{l + yf\ (156) 
and recall the trace evaluation of (|6^). We find 

ICl^i^ ^ pSm^n+l , (157) 

and the equation for a^"'' (t) then has the explicit form 

This differential recurrence relation is solved successively by 
4°)(t) = Trtel^l, 

.(atV)^(^-f)(e-^e--) (159) 

and so forth, where 



(ata) Uf\^\ =P , ... (160) 
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are the expectations values in the time-averaged steady state g'^^^'^ . 

Actually, Fig. | just shows this (a^a)^ for p = 0.7 and AT = 0.4, so that 
(a^a^^ = 1.75 and the approach to this asymptotic value is plotted for (^a'^a^^ = 

0, 0.35, 0.7 by the three dash-dotted curves. The time-averaged value of 

{a^a)^ is not well defined at t — 0, the instant of the first kick, any value in the 
range • • • 0.7 can be justified equally well. The memory of this arbitrary initial 
value is always lost quickly. 

4.2 Conditional and Unconditional Evolution 

Let us now be more realistic about the effect an atom has on the photon state. 
In fact, we have worked that out already in Sect, ^for the case of atoms incident 
in state — or in state ^: and a resonant Jaynes-Cummings coupling between 
the photons and the atoms. For the theoretical description of many one-atom 
maser experiments this is actually quite accurate, and it will surely do for the 
purpose of these lectures. 

Since the atoms should deposit energy into the resonator as efficiently as 
possible, they are prepared in the state. The net effect of a single atom is 
then available in (p5|), which we now present conveniently as 

Mgt^Agt+Bgt-gt^{A + B-l)gt (161) 

with 




(162) 



As the derivation in Sect. [l| shows, the term Agt corresponds to the atom emerg- 
ing in state and likewise Bgt refers to Accordingly, the respective proba- 
bilities for the final atom states ^: and — are 

prob(:=: ^ m) = Tr {Agt} , prob(:!: ^ :!:) = Tr {Bgt} , (163) 

and the probabilities p—, p^ that the state-selective detection of Fig. || finds 
the atom in ^: or are 

= Tr {Agt} , P^=V^ Tr {Bgt} , (164) 

respectively, where t]—,?]-^ are the detection efficiencies. 

The effect on the photon state of an atom traversing the cavity at time t can, 
therefore, be written as 

Qt ^ + V-Bgt + [(1 - v^)Agt + (1 - v-)Bgt] (165) 

where the three terms correspond to detecting the atom in state i^, detecting it 
in state and not detecting it at all. The probability for the latter case is 

prob(no click) = Tr {(1 - r/-)ylgt + (1 - V^)Bgt} = 1 - Tr {Cgt} , (166) 
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incident 
atom 



resonator 



atom detected 
in upper state 

atom detected 
in lower state 



(excited) 

Fig. 5. The final state of tfie exiting atom is detected: Is it or ^: ? 



where we recognize that 

Woh{:^ + prob(:!: ->:?:)= Tr {(^ + B)gt} = 1 (167) 

and introduce the click operator C, 

C ^r]-A + ri^B . (168) 

We take for granted that the atoms arrive with rate r at statisticaUy indepen- 
dent instants (Poissonian arrival statistics once more). The change of Qt brought 
about by a single undetected atom is 



Agt 



undetected atom 



{A + B-C)Qt 
1 - Tr {Cgt} 



Qt , 



(169) 



where the numerator is just the third term of (|165| ) and the denominator is its 
trace. We multiply this with the probability that there is an atom between t and 
t + dt, which is r d^, a nd with the probability that the at om escapes detection, 
which is given in (|166D and equal to the denominator in (169), and so get 



dt 



dgt 
dt 



undetected atoms 



rdt[iA + B-C)gt - gt+ Tt {Cgt} gt] 



We combine it with 1x2 



to arrive at 



dgt 
dt 



photon decay 



Cgt 



dt 



gt=[C + r{A + B- 1)] gt-r[C- Tr {Cgt}] gt 



(170) 



(171) 



(172) 



the master equation that applies between detection events. Owing to the term 
that involves the click probability Tt {Cgt}, this is a nonlinear master equation 
unless ?7— ~ rj-^ = rj when Tt {Cgt} = rj for all gt. Fortunately, the nonlinearity 
is of a very mild form, since we can write (172) as 



d_ 

dt 



gt = Cr^Qt - Tr {Cr^gt} gt 



(173) 



with the linear operator C^^ given by 



r{A + B -I) -rC = C + rM- rC 



(174) 
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and then solve it by 



Qt 



TrjeA^ft^o} 



(175) 



In effect, we can just ignore the second term in (173), evolve gt linearly from gt=o 
to e £it=Oj and then normalize this to unit trace. The normalization is necessary 
because £^ by itself docs not conserve the trace, except for r\— = 77-^ = 0. As 
we'll see in the fifth lecture, the normalizing deno minator of (175) has a simple 
and important physical significance; see Sect. 5.2. 

Note that there is a great difference between the situation in which the atoms 
are not observed (you don't listen) and the situation in which they are not 
detected (you listen but you don't hear anything). The evolution of the pho- 
ton field with unobserved atoms is the C = version of (172) that obtains for 
r\— = 1]-^ = 0, which is just the Scully-Lamb equation (145) with A4 of ( |l6l| ). 
It describes the unconditional evolution of the photon state. By contrast, if the 
atoms are under observation but escape detection, the nonlinear master equa- 
tion (172) or (173) applies. It describes the conditional evolution of the photon 
state, conditioned by the constraint that there are no detection events although 
detection is attempted. 

The difference between conditional and unconditional evolution is perhaps 
best illustrated in the extreme circumstance of perfectly efficient detectors, 77— = 
r]-^ = 1, when no atom escapes detection. Then (172) turns into (171) - as it 
should because "between detection events" is tantamount to "between atoms" 
if every atom is detected. More generally, if the detection efficiency is the same 
for and 77— = rj^ = ry, so that each atom is detected with probability 77, 
we have 

^^gt=[C+{l-v)rM]gt (176) 

for the evolution between detection events. This is the Scully-Lamb equation 
with the actual rate r replaced by the effective rate (1— 77)r, the rate of undetected 
atoms. 



4.3 Physical Significance of Statistical Operators 



Master equation (172) is the generalization of (176) that takes into account that 
ij: atoms are not detected with the same efficiency as atoms. This asymmetry 
may originate in actually different detection devices or - and this is in fact the 
more important situation - it is a consequence of the question wc are asking. 

Assume, for example, that the ^: detector clicked at t = and you want to 
know how probable it is that the next click of the i^: detector occurs between t 
and t + dt. Clicks of the detector are of no interest to you whatsoever. In an 
experiment you would just ignore them because they have no bearing on your 
question. The same deliberate ignorance enters the theoretical treatment: you 
employ (|172| ) with 77-^ = in the click operator ( |168D irrespective of the actual 
efficiency of the detector used in the experiment. Likewise if your question 



were about the next click of the detector you'd have to put 77— = in ( |168| ). 
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All of this is well in accord with the physical significance of the statistical 
operator Qt : it serves the sole purpose of enabling us to make correct predictions 
about measurement at time i, in particular about the probability that a certain 
outcome is obtained if a measurement is performed. Such probabilities are al- 
ways conditioned, they naturally depend on the constraints to which they are 
subjected. Therefore, it is quite possible that two persons have different statisti- 
cal operators for the same physical object because they take different conditions 
into account. 

Let us illustrate this point by a detection scheme that is simpler, and more 
immediately transparent, than the standard one-atom maser experiment speci- 



fied by A and B of (f 62). Instead we take 



AQt = ^ Qt ^ , 

i-(-ir'" i-(-ir'" 

l^Qt = ^ Qt ^ , (177) 

which can be realized by suitably prepared two-level atoms that have a non- 
resonant interaction with the photon field and a suitable manipulation prior to 
detecting ij: or What is measured in such an experiment is the value of 

(— 1)'*^'^, the parity of the photon state. Detecting the atom in state indicates 

even parity, (—1)"^'^ ~ 1, and a chck indicates odd parity, (—1)°"^" = —1. 

Now consider the four cases of Figs. ^ and |^. They refer to parity measure- 
ments on a one-atom maser that is not pumped (no resonant atoms are sent 
through) with v = 2 (an atypically large number of thermal photons for a mi- 
cromaser experiment) and r/A — 10. The plots show the period t = • • • 100/r 
of the simulated experiment. In Fig. ^ we see the parity expectation value as a 
function of t, and in Fig. ^ we have the expectation value of the photon number. 

The solid lines are the actual values, the vertical dashed lines guide 

the eye through state-reduction jumps, and the horizontal dash-dotted — • — • — 
lines indicate the steady state values 

^ i , (ata)^^^^-2. (178) 

On average, 100 atoms traverse the resonator in this time span, the actual num- 
ber is 108 here, of which 67 emerge in state ^: (even parity) and 41 in state 
— (odd parity). The final state is for 60% of the atoms on average, and 
for 40%. The values chosen for the detection efficiencies are r]— = 10% and 
77-^ = 15% so that each detector should register 6 atoms on average in a period 
of this duration. In fact, 7 zii clicks occurred (when rt=38.51, 44.80, 49.52, 53.07, 
72.05, 76.41, and 76.75) and 5 clicks (when rt=3.88, 85.81, 86.09, 94.12, and 
94.90). 

Experimenter Bob pays no attention to the even-parity clicks of the z^l de- 
tector, he is either not aware of them or has reasons to ignore them deliber- 
ately. Therefore he uses the nonlinear master equation ( |173| ) with 77— = and 
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77-^ — 0.15 for the evolution between two successive clicks of the — detector, and 
performs the state reduction 

whenever a — click happens. For example, to find the statistical operator ft=60/r 
and then the expectation values 

Bob, rt = 60: ((-1)"^") = 0.1920 , (0^0) = 1.787 , (180) 



he applies (|179| ) to the click at rt = 3.88, the only one between rt = 
and rt = 60. Bob's Qt=o is the steady state g^^^^ of the Scully-Lamb equation, 
consistent with his knowledge that the experi ment has been running long enough 
to have lost memory of its early history. For (177), g^'^^^ is actually the thermal 



state (|38|), here with 1^ — 2. Bob's detailed accounts are reported in Figs. ||(b) 
and g(b). 

Similarly, Figs, ^(c) and |^(c) show what Chuck has to say who pays no 
attention to clicks, but keeps a record of ^: clicks. He uses ( |173| ) with 77— = 0.1 
and 77-^ = for the evolution between two successive clicks and performs the 
state reduction 



for each ^: click. To establish 

Chuck, ri = 60 : ((-1)"^°^) = 0.3042 , (a^'a) = 1.599 (182) 



he has to do this for the four clicks prior to i = 60/r. Chuck uses the same pt^o as 
Bob, of course, because both have the same information about the preparation. 
And then there is Doris who pays full attention to all detector clicks. She 



uses (173) with rj— = 0.1 and 77-^ = 0.15 between successive clicks, does (181) 



for ij: clicks and (179) for clicks, and arrives at 



Doris, rt ^ 60 : ((-1)"^") = 0.2995 , (a^a) = 1.390 . (183) 

Her account is shown in Figs, ^(d) and 0(d). 

For the sake of completeness, we also have Figs, ^(a) and 0(a), where state 
reductions are performed for each atom, whether detected or not, and the 77— = 



77-^ = version of (173) - the "between atoms" equation ( |171| ) - applies between 
the reductions. What is obtained in this manner is of no consequence, however, 
because it incorporates data that are never actually available. 

Why, then, do we show Figs. |6|(a) and 0(a) at all? Because one might think 
that they report the "true state of affairs" so that 



aU atoms, = 60 : ((-1)"^") = 0.4818 , {a''a) = 1.189 (184) 



would be the "true expectation values" of (— 1)*^ and a^a a.t t — 60/r. And 
then one would conclude that the accounts given by Bob, Chuck, and Doris are 




Fig. 6. Different, but equally consistent, statistical predictions about the same physical system: Parity expectation value when parity 
^ measurements are performed on an unpumped resonator; see text 
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wrong in some sense. In fact, all three give correct, though differing accounts, 



and the various predictions for t — 60/r in ( |180| ), ( |182| ), and (183) are all sta- 
tistically correct. For, if you repeat the experiment very often you'll find that 
Bob's expectation values are confirmed by the data, and so are Chuck's, and so 
are Doris's. 

But, of course, when extracting (a^a)t^go/r' ^^^^ from the data of the very 
many runs, you must take different sube nsem bles for checking Bob's predictions. 



or Chuck's, or Doris's. Bob's prediction (180) refers to the subensemble charac 



terized by a single click at rt = 3.88 and no other click between rt — and 



rt = 60, but any number of z^l clicks. Likewise, Chuck's prediction (182) is about 
the subensemble that has m clicks at rt=38.51, 44.80, 49.52, and 53.07, no other 
li: clicks and any number of clicks. And Doris's subensemble is specified by 
having this one click, these four clicks, and no other clicks of either kind. 
By contrast, no experimentally identifiable ensemble corresponds to Figs, ^(a) 
and 0(a) ; they represent sheer imagination, not phenomenological reality. 

In summary, although it is the same physical object (the privileged mode of 
the resonator) that Bob, Chuck, and Doris make predictions about, they regard 
it as a representative of different ensembles, which are respectively characterized 
by the information taken into account. This illustrates the basic fact that a 
statistical operator is just an encoding of what we know about the system. 
Depending on the question we are asking, we may even have to ignore some 
of the information deliberately. The appropriate Qt is the one that pays due 
attention to the pertinent conditions under which we wish to make statistical 
predictions. Different conditions simply require different statistical operators. 
That there are several git's for the same physical object is then not bewildering, 
but rather expected. 

The ensembles that Bob, Chuck, and Doris make statements about need not 
be, and usually are not, real ensembles created by many repeated runs of the 
experiment (or, since the system is ergodic |3^ , |3^ , perhaps a single run of very 
long duration). Rather, they are Gibbs ensembles in the standard meaning of 
statistical mechanics: imagined ensembles characterized by the respective con- 
straints and consistent weights. The constraints are duly taken into account by 



the nonlinear master equation (173) 



We should also mention that Bob's Qt is consistent with Chuck's and Doris's 
by construction. If we didn't know how they arrive at their respective statisti- 
cal operators we might wonder how we could verify that the three gt's do not 
contradict each other. Consult p^,pq| if you find the question interesting. 



Homework Assignments 

17 Show that 



y y^gW = -1±^ ( 1 + y ) (185) 
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is the V > generalization of (|5^), then use this and (jsj) to find the ly > 
version of (157). 



18 State the > version of (158) and solve the equations for n = 0, 1, 2. 

19 Take the Scully-Lamb limit of (156), that is: T after putting p — rT, 

and find the steady state values of all al^'' 's. 

oo 

20 Evaluate p^^"^ = ^ "1"'^!°^ for these d'h^'s. You should obtain 



Tl = 



1 



1 



a'^a \ i/AJ 



-r/A 



(186) 



In which sense is this yet another generating function for the 's? What 
do you get for i' 01 

21 Verify ( p7| ) for B of ( |l62| ). 

22 Show that 

Mf{a^a) = sin(^0\/aF^^ ^/(a^a-1)- sin (^^ata + l^ V(a^a) (187) 

for M of (161), then combine this with ( ^5| ) to find the steady state of the 
Scully-Lamb equation (145). You should get 



1 



r/A 



(188) 



as the steady state of the one-atom maser. What is the physical significance 
of the normalizatio n fa ctor iV? How do you determine it? 

23 Differentiate gt of (175) to verify that it solves (173). 

24 In addition to Bob, Chuck, and Doris, there is also Alice who pays attention 
t o al l detector clicks but doesn't care which detector fires. Which version of 
( |l73| ) does she employ, and how does she go about state reduction? 

25 (a) Since you found Sect. 4.3 very puzzling, read it again. 

(b) If you are still puzzled, repeat (a), else proceed with (c). 

(c) Convince your favorite skeptical colleague that there can be different, 
but equally consistent, statistical operators for the same object. 



5 Fifth Lecture: Statistics of Detected Atoms 

A one-atom maser is operated under steady-state conditions. What is the prob- 
ability to detect a atom between t and t + dtl It is 

rdtr]-TT ^^Ag^'"''^ , (189) 

the product of the probability rdt of having an atom in this time interval (r > 
is taken for granted from here on) and the probability rj— Tr {Ag^^'^^ } that the 
detector clicks if there is an atom. Similarly, 

rdtT]^Ti ^^Bg^''^^ (190) 
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is the probability for a — click between t and t + dt. What multiplies rdt here 
are the traces of the first and the second term of (165), respectively. The third 
term gave us the no-click probability ( |166|). 

The a priori probabilities (189) and ( |19C| ) make no reference to other detec- 
tion events. But the detector clicks are not statistically independent. For the 
example of the parity measurements of Sect. 4.3, you can see this very clearly in 
Fig. ^(d), where the even-parity clicks and the odd-parity clicks come in bunches. 
This bunching is easily understood: a zii click is accompanied by the state re- 
duction (|l8l| ) so that Qt is an even-parity state immediately after a ^: click and, 
therefore, the next atom is much more likely to encounter even parity than odd 
parity. The detection of the first atom conditions the probabilities for the second 
atom, and for all subsequent ones as well. 



5.1 Correlation EYinctions 



Perhaps the simplest question one can ask in this context is: Given that a 
atom was detected at i = 0, what is now the probability to detect a — atom 
between t and t -\- di? Sinc e de tection events that occur earlier than t are not 
relevant here, we must use (173) with rj— = rj-^ = to propagate pt=o to Qt, 

Qt = e^o*et=o (191) 



where 

ro-A = C + r{A + B -I) (192) 

r]——r]^—Q 

is the Liouville operator of the Scully-Lamb equation that determines the steady 
state, 

£0^^^''''^ = , lim e^'^^Qt=o = q'^''''^ • (193) 



The denominator of (175) equals unity here because Co conserves the trace, 

Tr { e^«*^)t=o} = Tr {Qt^o} = 1 . (194) 

With 



Tr 

accounting for the initial click at t — (which, under steady-state conditions, 
is really any instant), the probability for a click at t • • • t + dt is then 



rdtrj-^ Tr ^Bc~"''Qt^ 



.Cot, 



= rdt 



Tr \b e'^ot j^q(ss)^ 
Tr 



(196) 



If t is very late, the memory of the initial state gets lost and ( |196| ) becomes equal 
to the a priori probability (190). The ratio of the two probabilities, 



Tr|Se^o*y^^?(^^)} 
Tr {Bq^^^)] Tr 



(197) 
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Fig. 8. Correlation functions for parity measurements on an unpumped resonator. The 

solid line shows the even-even correlation function G——(t), the dashed 

line the odd-odd correlation function G^^{t), and the dash-dotted — • — ■ — hue the 
cro ss co rrelati ons G— — (t) — G-^—{t). The plot is ior u — 2 in the expressions given 



in (204) and (205) 



is the correlation function for clicks after ij: clicks. There are no correlations if 
G— ^ = 1, positive correlations if G— ^ > 1, negative correlations if G— ^ < 1. 
The bunching of Fig. |^(d) shows strong negative zii— correlations (odd after 
even) at short times. 

Likewise, if we ask about z^i clicks after clicks, we get 



Tr \^Ae^otQg{ss))^ 
Tr Tr {Bg<-^^)} 



(198) 



and we have 



G^^it) = 



Tr ^^Ae^otJ^g(ss)^^ 
^TY{^4^M}]^ 



G^^{t) = 



Trjee^o^e^"'')} 
[Tr{6eM}]^ 



(199) 



for the correlation functions of clicks of the same kind. Note that the detection 



efficiencies rj—^ri-^ do not enter the correlation functions (197)-(199). This is a 
result of normalizing the conditional probabilities to the unconditional ones. But 
see homework assignment 

As an example we take once more the parity measurements of Figs. ^ and ^ 
for which A, B are given in ( p77| ) so that Co=C. Then is the thermal state 
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of ( |38| ) and the a priori probabilities are 

1 1 



Tr{^^?M} = i + i((-l)«'-) 
Tr{Se(-)} =i-i((-ir'") 



(ss) 



V 

2iy + l ■ 



(200) 



For 1/ = 0, there are no odd-parity chcks of the detector and, therefore, we'll 
only consider the case of > 0. Since 



2i^ + l V 



we have 



and 



(201) 



(202) 



(203) 



(197)-(199) are then easily evaluated and we obtain 



where X{t) is given in ( |49| ) with A(0) = (2;/ + l)/(j^ + !)■ The numerators of 

(204) 



G^-{t) = G^-{t) = 1 - l + (2z/+l)2(e^* 



1) 



and 



G--{t) = 1 
G^^it) = 1 



1/ 



V 



l + (2j/+l)2(c^*-l) 



l + (2z/+l)2(e^*-l) 



(205) 



The cross correlations of ( 204 ) vanish at i = (no even-parity click immediately 
after an odd-parity click and vice versa) and increase monotonically toward 1. 
The two same-click correlation functions of (205) are always larger than 1 and 
decrease monotonically from their i = values 



G'--(t = 0) 



Z/+ 1 



2z/+ 1 



(206) 



which exceed unity and so confirm the bunching observed in Fig. ^(d). For 
1/ = 2, the parameter value of Figs. ^and|^, the correlation functions are plotted 
in Fig. I 

Data of actual measurements of correlation functions for atoms emerging 
from a real- life micromaser are reported in ||3^ , for example. Theoretical values 
for some related quantities, such as the mean number of successive detector clicks 
of the same kind, agree very well with the experimental findings. 
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5.2 Waiting Time Statistics 

Here is a different question: A click happened at f = 0, what is the probability 
that the next ^: clicks occ urs between t and t + dt? In marked contrast to the 
question asked before (191), we are now not interested in any later click, but in 
the next click, and this just says that there are no other ^: clicks before t. Since 
we ignore deliberately all clicks at intermediate times, we have to use (173) 
with T]-^ = in the click operator C of (168). 
We introduce the following quantities: 

rdtp—{t) = probability for a ^: click at t ■ ■ - t + dt , 

Pno—{t) " probability for no click before t , 

dt Pnoxt— (i) — probability for the next click to happen ai t ■ ■ -t + dt . 

(207) 

Since 1 — rdtp—{t) is then the probability that there is no i^: click between t 
and t + dt, we have 



Pno- {t + dt) = pno- (t) [I - rdt p- (t)] 

, Id 

or PizW = --T7lnpno-(i) , 
r dt 



(208) 



and 



Pno-it + dt) ^ Pno-{t) ~ dt P„cxt-(i) 

or Pnc^t-it) = --jrPno-it) 
dt 



(209) 



is another immediate consequence of the significance given to these quantities. 
We knowp— (t) from (166) and (175), 



P^it) 



Tr|c e^'J^et^o} 
Tr{e^'7^e(=o} 



(210) 



where 



C = 



and 



gt=o = 



Tr{Sg('^^)} 



(211) 



in the present context. As required by (208), the right-hand side of (210) must 
be a logarithmic derivative and, indeed, it is because the identity 



^ Tr { e^^^gt=o} = Tr {c„ e^^^Qt=o} 



Tr 



[{Co-rC)e^^^^gt=o} 
-rTr{Cc^"*£)i=o} , 



(212) 
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0.005 



Fig. 9. Waiting time distributions for parity measurements on an unpumped resonator. 

The solid lines show Pncxtit)/r for clicks after z^i clicks (left) and clicks 

after clicks (right) as functions of rt. The dashed lines refer to clicks after 

clicks (left) and i^: clicks after clicks (right). In these logarithmic plots, straight 
lines would correspond to the Poissonian statistics of uncorrelated clicks 



which uses the trace-conserving property of Cq , implies 

1 d 



It follows that 



Pno-{t) = Tr{e^"*gt=o} 



and then (^09|), (|212|) give 

P„oxt-(t) ^rTrjCe^^^ft^o} 



(213) 
(214) 

(215) 



The "important physical significance" of the denominator in ( |173| ) that was left 
in Umbo in Sect. 4.2 is finally revealed in ( ^14[ ): it is the probability that no 
atom is detected before t. Since an atom is surely detected if we just wait long 
enough, the limit 

Pno—{t) — > as t — > oo (216) 

is a necessary property of Pno—it). As a consequence, all eigenvalues of Cfj must 
have a negative real part. 

Putting all things together we obtain 



-Pnext-(0 = rri— 



Tr{^e(^o-'^'?--^)*6e(^^)} 



Tr {SpM} 



(217) 



for the waiting time distribution for the next click after a click. Analogous 
expressions apply for the next — click after a ^: click, the next z^i click after a 
^: click, and so forth. As a basic check of consistency we consider the situation 
in which A and B arc just multiples of the identity. 



Agt = qgt , Bgt = (1 - q)gt with < g < 1 . (218) 
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Then the detector clicks are not correlated at all and the waiting time distribu- 



tion (217) should be Poissonian, 

P„cxt-(i) = rr^-qe-'^'^^l^ , (219) 



and this is indeed what we get from ( piT] ) for ( |218| ). 

Figure ^ shows the waiting time distributions to the parity measurements of 
Figs. Other examples are presented in some figures of |3^ . 

5.3 Counting Statistics 

Yet another question is this: What is the probability Wn{t) for detecting n atoms 
in state ^: during a period of duration t? We pay no attention to :Tz clicks and, 



therefore, have 77-^ = in the nonlinear master equation (173) for the evolution 
between clicks. 



Probability wo(i) is the no-click probability of Sect. 5.2, 



Wo = Tr I e'^"* | with C = rj-A , (220) 



where gt=o — Q^^^'^ in ( |214| ) is appropriate now. The one-click probability wi{t) 
is given by 



(221) 

where the last factor is the probability wo{t') for no click before t', the first factor 
is the probability for no click after t', and 

rdt' TT{Cgt=t'-o} (222) 



is the probability for a click at t' . In accordance with (175) and (181), the 
statistical operators just before and after the click at t' are 

Tr I e^"* e^*^^) > 

and 

_ Agt=t'-o _ Cgt=t'-Q , . 

gt=t'+Q — rr. r A 7 — 7 ' l^^^j 

Tt {Agt=t'-o} Tt {Cgt=t'-o} 

respectively, so that 

gt=t'+oTT{Cgt=t'-o} Tr { e^''*'^.^''^) } = C e^'-^'e^'''') , (225) 
and a remarkable simplification happens, inasmuch as 

wi{t) = J dt' Tr I e^''^* ~ *')rCe^''*'g>("^'} (226) 
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involves but a single trace as the equivalent replacement of the product of three 
traces with which we started in (^2l|). 

Upon writing (226) as a double integral 

wi{t) = J dti J dtoS{to + h-t)TT \^e^^^^rCe^^^Og'^''''^j (227) 

it is reasonably obvious (and can be demonstrated by a simple induction) that 



Jo Jo Jo 

Tr I e^^^'VCe^'^'^^-irC • • • e'^'f^irC e^^^V'""' } (228) 



or 



Wn{t) = Tr [Wn{t)g(''^} . (229) 
The operator Wn{i) thus introduced, 



oc 



Wn{t) = I dtn ■■■ I dto S{to + --- + tn-t) e^'i^^rC e^"*"-! • • • rC e^"*" , 

(230) 







obeys the recurrence relation 

Wnit)^ f dt'Wo{t-t')rCWn-i{t') for n = 1,2,3,... (231) 
that commences with 

Wo{t) = e^"'^ . (232) 
As always, we'll find it convenient to use a generating function, 

oo 

n;(x,t) = ^x" W„(i) . (233) 

n=0 

The recurrence relation for the yV„ {t) 's then turns into an integral equation for 
W(x,t), 

W{x,t)^Wo{t)+ [ dt'Woit~t')xrCW{x,t') . (234) 
Jo 

The standard technique for handling such finite-range convolutions utilizes 
Laplace transforms because the familiar identity 

dte~"<* dt'f{t-t')g{t')^ dte-^^f{t) dt' e''^^ g{t') (235) 
Jo Jo Jo 

leads to a factorization. With 

dt e~^^Wo{t) = (7 - Crj)'^ (236) 
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this gets us to 

poo 

/ dt e-^^W{x,t) = [1 - ~ £,,y\rC]'\^ - 
Jo 

= (7 - Cr, - xrC) , 
and the inverse Laplace transform is elementary, 



(237) 



As we noted at ( pi6D , all eigenvalues of £^ have negative real parts, and so 
the Laplace transform ( ^36[) of Woit) surely exists for 7 > 0. The same remark 
applies to yV(x, t) for —(1 — Tj—)/ri— < x <1 because 



+ xrC = £0 — (1 — x)rC = Cq — r{l — x)r]—A 



(239) 



is with Tj— replaced by (1 — x)rj— so that Cr, + xrC is just another operator 
of the same kind as £^ if the "effective detection efficiency" (1 — x)ri— is in the 
range • • • 1, which restricts x to the range stated above. 
After putting things together we obtain 



^ x'^Wnit) = Tr { e(^'' + ^^C)i^(ss)| 



(240) 



as the generating function for the counting probabilities Wn (t) ■ In view of what 
we observed above about £^ + xrC, the right-hand side of ( |24Cl| ) is equal to 
the no-click probability WQ{t) for detection efhciency (1 — x)rj—. Accordingly, 
Wo(t) determines all w„(t) through its dependence on 77—. As an immediate 
consequence of this observation, we get a statement about the moments of the 
counting statistics. 



k 00 



n=0 ^'^^ ^ ri=0 

For fc = 0, we check the normalization, 



ri_->0 



^u;„(t) = Tr{e^o*^,(-)} = l; 

n=0 

for fc = 1, we get the average number of ^: clicks, 

Vn«;„(i)= / dt' Tr|e^"(*-*')r?7-yle^"*'f?("^)| 
n=o Jo ^ ^ J 



(241) 



(242) 



(243) 
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which can be understood as a statement about the ergodicity of the process |Q ; 
and for fc = 2, we learn something about the variance of the counting statistics, 

oo 

n{n - l)w„{t) 

2 J dt' J dt" Tr I e^«(* " ^'^rij-Ae^''^^' ^ ^"^r-q-Ae^"^" g^''^^ 
{rj-rtf Tt {^AE{Cot)Ag^''^} , (244) 



n=0 



with 



E{y) = -{ey-l-y) 



(245) 



In traces of integrals such as ( |243D and ( p44| ), the exponential on the left and 
on the right can be ignored because Cq is trace conserving and g'^^^^ is its right 
eigenvector to eigenvalue zero. 



Here, too, we get Poissonian statistics for ( 218 ), namely 

^ — ^ rj ' 



n=0 

for which 



n=0 



{rrj—qt)^ 
k\ 



In particular, we note that 

oo 

n{n - l)wn{t) = 



n=0 



H 2 



nWn{t) 

.n=0 



(246) 



(247) 



(248) 



holds for the Poissonian counting statistics ( p46 

A convenient, yet rough, measure for the deviation from Poissonian statistics 
is the so-called Fano-Mandel factor Q{t), 



Q = 



n=0 



ri=0 



(249) 



n=0 



a normalized variance. The normalization is such that Q = Q for Poissonian 
counting statistics, as one verifies easily with ( |24^ ). For — 1 < Q < one speaks 
of sub-Poissonian statistics, and of super-Poissonian statistics for Q > Q. 
For the count of m clicks, (|243D and (|244|) give 



Ti {AE{Cat)Ag'^''^^] 
Tr {y^pM} 



Tr 



[Ag^^^^} 



(250) 
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If 



0.12 



0.09 



0.06 



0.03 




Fig. 10. Fano-Mandel factor for the counting statistics in pari ty measurements on an 
unpumped resonator. The solid hnc s hows Q—{t) of (252), the horizontal dashed 



line is the asymptotic value of ( 253 ) 



We use the damping bases to get a tractable numerical expression, 

n— A:— — oo 



Tr 



^^^^^ 

where the n = 0,k = term is removed since £»q°'' = g'''^^\ ^g"'' — 1, Aq^"* = 0, 
and E{y = 0) — 1. For the parity measurements of Figs. ||-^ one can evaluate 
the sum and gets 



rj—r 



dt' 



■In 



1 



{iy + l){2iy+l) Jo 2At 
which is plotted in Fig. ^ For At <^ 1 and At ^ 1 the limiting forms 



2v+l 

\tl{2v + 1) 



Tj—rt for At < 1 , 

77— r/A for At > 1 , 



(252) 



(253) 



{v + l){2v+l) 



obtain, as is confirmed by Fig. |lO|. Other examples of Fano-Mandel factors for 
counting statistics are presented in some figures of [^5| . 

The pioneering measurement in 1990 of atom statistics in a real-life micro- 
maser experiment is reported in ]40[ | and linked to the photon counting statistics 
in H] . Measured Fano-Mandel factors from some later experiments can be found 
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Homework Assignments 

26 What is the correlation function for cUcks of either kind, that is: without 
caring if it's a chck or a chck? 

27 The two c ross- correlation functions are identical in the example of ( ^OOD - 
(206), see ( 204 ). Is this always the case? 



28 Since the next click is bound to come sooner or later, consistency requires 
that Pncxtjt ) is normalized to unit integral. Show that this is indeed so for 

P„ext(0 of (111). 

29 A ^: click happens at t = 0. What is the probability that the next click is a 

click? 

30 Use the methods of Sect. to find an expression for the average waiting 
time between successive m clicks, between successive clicks. 

31 Determine the short-time behavior of the various Pnext's of Fig. ^ and com- 
pare with the plots. 

32 Show that the exponential function e"^ of an operator F responds to varia- 
tions SF in accordance with 

S = [\t e^^SF e^^ " , (254) 
Jo 

which epitomizes all of perturbation theory. 



33 Use ( P54[ ) to extract Wi{t) and W2(i) from ( |24C| ). Compare with ( p28[ ). 

34 Consider the probability Wnmit) of detecting n atoms in state z^i and m 
atoms in state — during a period of duration t. Show that 



E 

n.m— 



x"y'^Wnm{t) = Tr{e(^" + ^'^''-'4 + ^^'^''-^)*^^(''^)} (255) 



is the appropriate generalization of the generating function (240) 



35 For A of (177), evaluate the k = traces in (U5w and so confirm (|252|) 



36 Find the leading correction to the approximation given in (|253|) for At^l. 
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Appendix 

Here are some facts about special functions that are useful for homework as- 

signments ^ and ^. The expansion of ( [73| ) in powers of e"^"^* cx [Q;(t)/a*(i)] ^ is 
done with the aid of 

oo 

^l^x{y-l/y)^ J2 y^J.ix), (256) 
k— — oo 

the most important generating function for Bessel functions of integer order, 

(-1)™ ^l^^^|fe|+2m 

(257) 



j,{x) = (-i)^-j_,(x) = ii'=i-'=j|.|(x) = ii'^i-'^ , ii^hi (H' 

' ' ^-^ m\ [m + \k\)\ ' 



m=0 



They in turn act as a generating function for Laguerre polynomials, 

oo „ 

^^2^) = (.,)^|fcle-^^^-^Lf l)(.) . (258) 

n— ^ II/" 

After a suitable Laplace transform this becomes 

oo 

(1 + y)-\^\-^ e^ - ^(~y)"L(l'=l)(a;) , (259) 
which is another useful generating function for the Laguerre polynomials 

\m + fc / m! 

m— ^ 

The integral relations 



1 1 /"^ 1 

(Ifc|)(a;)^ e^x-2W / d2/e-V+^l^'lj|fe|(2V^) (261) 



and 



/ dye-2^J|,.|(2ViI^)j|fc|(2V^y^) = e-(^ + ")lfe(2V^) (262) 
Jo 
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are worth knowing, where 

Ife(x) =r%(ia;) (263) 

are modified Bessel functions of integer order. As a preparation for homework 
assignment ^ you might want to derive first 

(264) 

by combining some of these relations fittingly. 
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